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The approximate renormalized stress-energy tensor of the quantized mas- 
sive conformally coupled scalar field in the spacetime of electrically charged 
nonlinear black hole is constructed. It is achieved by functional differentiation 
■ of the lowest order of the DeWitt-Schwinger effective action involving coinci- 

<^> . dence limit of the Hadamard-Minakshisundaram-DeWitt-Seely coefficient 03. 

The result is compared with the analogous results derived for the Reissner- 
Nordstrom black hole. It it shown that the most important differences occur 
' in the vicinity of the event horizon of the black hole near the extremality 

\& ■ limit. The structure of the nonlinear black hole is briefly studied by means of 

. the Lambert functions. 
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O I. INTRODUCTION 

o ■ 

For the quantized massive fields in the large mass limit, i. e., when the Compton length 
is much smaller than the characteristic radius of a curvature, the nonlocal contribution to 
the effective action can be neglected, and the series expansion in m~ 2 of the renormal- 
>• ■ ized effective action, Wr, may be easily constructed with the aid of the DeWitt-Schwinger 
method [1,2]. As the renormalization prescription requires absorbtion of the first three terms 
of the series into the quadratic classical gravitational action, the n — th term of Wr is propor- 
tional to the integrated coincidence limit of the Hadamard-Minakshisundaram-DeWitt-Seely 
coefficient (HMDS) a n+3 . Unfortunately, the complexity of the coefficients [a n ] rapidly in- 
creases with n, and consequently one expects that the applicability of the series expansion 
is confined to a first, perhaps a second nonvanishing term. 

Having constructed a first order Wr, the renormalized stress-energy tensor (which is 
the most important characteristics of the quantized field in the curved spacetime) may be 
obtained in a standard way, i. e., by functional differentiating the constructed effective action 
with respect to the metric tensor. This method has been successfully applied in calculations 
of the approximate renormalized stress-energy tensor of the quantized massive scalar, spinor, 
and vector fields in the vacuum type-D geometries by Frolov and Zel'nikov [3-7]. 

A different method, based on the WKB approximation of the solutions of the massive 
scalar field equation in a general spherically-symmetric spacetime, and summation thus ob- 
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tained modes by means of the Abel-Plana formula, has been invented by Anderson, Hiscock, 
and Samuel and applied in the context of the RN spacetime [8]. Their method is equivalent 
to the Schwinger - DeWitt expansion: to obtain the lowest (i. e. m~ 2 ) terms, one has to 
use sixth-order WKB approximation. Moreover, numerical calculations reported in Ref. [8] 
confirmed that the DeWitt-Schwinger method provide a good approximation of the renor- 
malized stress-energy tensor of the massive scalar field with arbitrary curvature coupling as 
long as the mass of the field remains sufficiently large. This approach and its modifications 
has been employed in various contexts in Refs.[9-15]. 

Recently, extending the results of Frolov-Zel'nikov, we have constructed a general formula 
describing the approximate renormalized (T^) ren of the quantized massive scalar, spinor, 
and vector fields in arbitrary spacetime. The results have been presented for the class of 
geometries with vanishing curvature scalar, and subsequently applied in the spacetime of the 
RN black hole and in the spacetime that could be obtained by expanding its near horizon 
geometry into a whole manifold [16]. Our formulas allow, in principle, to determine the 
renormalized stress-energy tensor of the massive field once the line element has been chosen, 
although the specific calculations may be very tedious. For the quantized massive scalar 
field with arbitrary curvature coupling in the RN spacetime we have reproduced the results 
of Anderson, Hiscock, and Samuel; neutral spinor and vector fields have not been discussed 
earlier. 

In this paper we shall extend the analyses of Ref. [16] to the general geometry and con- 
struct the renormalized stress-energy tensor of the massive quantized scalar field obeying 
the equation 

(-□ + £R + m 2 ) = 0, (1) 

where £ is the coupling constant and m is the mass of the field. Since the background 
geometry is general, the most direct approach is to use the first nonvanishing term of the 
renormalized effective action. The advantage of this approach lies in the purely geometric 
nature of the approximation that reflects its local nature. Although the constructed result 
is rather complex, we shall present it in its full length, because it provides the generic for- 
mula from which the renormalized stress-energy tensor in some physically interesting cases 
may be easily obtained. As the effective action of the quantized massive scalar field differs 
form the analogous actions constructed for fields of higher spins only by numerical coeffi- 
cients, one can generalize presented results to fields of other spins. It should be emphasized 
however, that the method has obvious limitations, and, when applied to rapidly varying or 
strong gravitational fields it breaks down. Moreover, its massless limit is contaminated by 
nonphysical divergences. 

Our general formulas will be employed in the calculations of (T^) ren in the geometry 
of the electrically charged black hole, being an exact solution of the coupled system of the 
Einstein equations and the equations of the nonlinear electrodynamics recently proposed 



by Ayon-Beato and Garcia (ABG) in Ref. [17], to which the reader is referred for physical 



motivations and technical details. Their exact solution is characterized by the electric charge, 
e and the mass M, and may be elegantly expressed in terms of the hyperbolic functions. 
An important and interesting feature of this solution is its regularity as radial coordinate 
tends to zero. We shall show that the structure of horizons of the ABG solution may by 



studied by means of the Lambert function [18], allowing analytical treatment of the vacuum 
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polarization effects on the event horizon. At large distances their solution behaves as the 
RN solution. For small and intermediate values of the ratio \e\/M, the location of the event 
horizon, r + , is close to the location of the event horizon of the RN black hole; significant 
differences occur near the extremality limit. It would be, therefore, interesting to analyze 
how the similarities of the line elements are reflected in the behavior of the renormalized 
stress-energy tensors. 

The renormalized effective action of the massive scalar field involves the terms that are 
proportional to the first and third power of £ — 1/6. As the curvature scalar of the RN 
spacetime vanishes, (T^) ren of the massive scalar field naturally divides into the part that 
describes pure conformal coupling and an additional local part that is multiplied by a factor 
£ — 1/6. On the other hand however, the curvature scalar of the ABG geomerty does not 
vanish, and the structure of the effective action indicates that the renormalized stress-energy 
tensor of the massive scalar depends on the constant £ in a more complicated way. Since 
the conformal coupling leads to massive simplifications, one expects that the similarities in 
the renormalized {T^) ren (if any) would appear mainly in this case. 

II. THE RENORMALIZED STRESS-ENERGY TENSOR OF THE QUANTIZED 

MASSIVE SCALAR FIELD 

The renormalized effective action constructed for the quantized scalar field satisfying 
equation (1) is given by 

Wren = TTTj \ 2 / d ^ £ T^Kl (2) 

where [a n ] is the coincidence limit of the n— th HDSM coefficient. The first three coefficients 
of the DeWitt-Schwinger expansion, ciq, ai, and 02, which contribute to the divergent part 
of the action have to be absorbed in the classical gravitational action by renormalization of 
the bare coupling constant. 

As the rigorous asymptotic analysis of the fundamental solution is restricted, in general, 
to the heat operator of a parabolic type, we tacitly assume that all steps that are necessary 
in construction of the first-order renormalized stress-energy tensor have been carried out in 
an analytically continued Euclidean spacetime. The analytic continuation to the physical 
space is performed at the last stage of the calculations. 

Calculation of the HDSM coefficients is an extremely laborious task, and their exact form 
for n > 5 is unknown. The coefficient [a 2 ], which is proportional to the trace anomaly of 
the renormalized stress-energy tensor of the quantized, massless, and conformally invariant 
fields, has been calculated by DeWitt |lj. The coincidence limit of the coefficient 03 has 
been obtained by Gilkey [19,20] whereas the coefficient [04] has been calculated by Avramidi 
[21-24]. 

Restricting ourselves to the terms proportional to m~ 2 , one has 

Wren = ^>n~~2 2 / ^V^M, (3) 

with 
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where 



M = % + £5, (4) 
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Since the coincidence limit of the coefficient 0,4 is much more complex one expects that 
using it in the calculations of the approximate renormalized stress-energy would be a real 
challenge. However, it still could by of use in the simpler analyses of the field fluctuation, 
(0 2 )ren- Substituting (4-6) into (3), integrating by parts and making use of elementary 
properties of the Riemann tensor, one can reduce the number of terms in the renormalized 
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where en, are numerical coefficients that stand in front of the geometrical terms. 
The renormalized stress-energy tensor is given by the standard relation 

2 " W£J = (T ab ) ren . (8) 
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Functionally differentiating the renormalized effective action with respect to the metric ten- 
sor, performing simplifications and necessary symmetrizations, after rather long calculations, 
one has 
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As there are numerous identities involving the Riemann tensor, its covariant derivatives 
and contractions, the form of {T b ) ren is, of course, not unique and depends on adopted 
simplification strategies. Here we presented our results in the form that we have found 
useful in the further calculations. It should be noted that the resulting renormalized stress- 
tensor of the massive scalar field depends on the coupling constant in a complicated way, 
and in a general spacetime it divides naturally into four terms 

(T b a )ren = ^ (t ~ \) ^ + \ (t ~ ^ ^ + (* " ^ ^ + T^, (19) 



where 



Til)ab = ^- 1/2 iL iW5 - w ^ (20) 



T(2)ab = w«-*-L w - (2i) 



T 0W> = _±_ g -W* w ^ (22) 
96vr 2 m 2 5g ab 



T^ ab = ^r~ 2 g- 1/2 (fiJ-Wi + a*w 2 + , (23) 

96ir 2 m 2 \ bg ab bg ab ^ dg ab J 



and 



e = M-£>- (24) 

Inspection of eqs. (9-18) shows that variational derivatives of Wi and W 3: with respect to 
the metric tensor vanish in R = geometries, and, additionally, that of W 2 , W4, W e , and 
W7 vanish for the Ricci-flat geometries. Moreover, one has important simplifications of the 
general stress-energy tensor for the conformally coupled massive fields as there is no need 
to compute T^ ab , T^ a \ and T^ ab . Finally we observe that the analogous expression of 
the stres-energy tensor of the quantized massive spinor and vector fields differs only by the 
numerical coefficients a*. Inserting appropriate coefficients listed in the Table I into (7), 
one may easily generalize our discussion to the fields of higher spins. Note however, that to 
obtain the appropriate result for the neutral spinor field one has to multiply the renormalized 
effective action by the factor 1/2. 



III. ELECTRICALLY CHARGED NONLINEAR BLACK HOLE 

As is well known the Reissner-Nordstrom line element is the only static and asymptot- 
ically flat solution of the Einstein- Maxwell equations representing a black hole of mass M 
and electric charge e. The appropriate line element has the form 
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ds 2 = -U(r)dt 2 + V~ l {r)dr 2 + r 2 (sin 2 6d(j) 2 + d6 2 ) , (25) 

where the metric functions U(r) and V(r) are given by 

2M p 2 

U(r) = V(r) = 1 + -. (26) 

Because of its simplicity the RN solution may be studied analytically; for e 2 < M 2 the 
equation <7 o = has two positive roots 

r± = M ± (M 2 - e 2 ) 1/2 , (27) 

and the larger root represents the location of the event horizon, while r_ is the inner horizon. 
In the limit e 2 = M 2 horizons merge at r = M, and the RN solution degenerates to the 
extremal one. The singularity of the RN line element that one encounters at r = is a 
non-removable curvature singularity, while those at r± are merely spurious singularities that 
may be easily removed by a suitable choice of coordinates. 

Recent interest in the nonlinear electrodynamics is partially motivated, beside a natural 
curiosity, by the fact that the theories of this type frequently arise in modern theoretical 
physics. For example they appear as effective theories of string/M-theory. Moreover, one ex- 
pects that it should be possible to construct solutions to the coupled system of the Einstein 
field and equations of the nonlinear electrodynamics, which may be interpreted as repre- 
senting globally regular black hole geometries, avoiding thus the singularity problem. As 
the nonlinear electrodynamics in the weak field limit coincides with the Maxwell theory, one 
expects that the appropriate solution should approach at large distances the RN solution. 

An interesting solution of this type, representing spacetime of the regular black hole with 



mass M and charge e has been constructed recently by Ayon-Beato and Garcia fl? |. The 
appropriate line element is given by (25) with 



U(r) = V(r) = 1 - — 



, ( e ' \ 
tanh — — — 

\2Mr ) 



(28) 



For e = the ABG solution reduces to the Schwarzschild solution; for small values of the 
charge it differs from the Reisner- Nordstrom solution by terms of order 0(e 6 ). At large dis- 
tances the metric structure of (28) also closely resembles that of the RN solution. Indeed, 
expanding U (r) in a power series one concludes that the ABG solution behaves asymptoti- 
cally as 

2 /If p 2 p 6 1 

U(r) = V(r) = !-_ + _- + O(-). (29) 

Instead of referring to numerical calculations at this stage of analyses of the ABG geom- 
etry, we show that although the metric coefficient U(r) is a complicated function of r, the 
location of the horizons may be elegantly expressed in terms of the Lambert functions 



Indeed, making use of the substitution r = Mx and e 2 = q 2 M 2 , and subsequently introduc- 
ing a new unknown function W by means of the relation 

4o 2 
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one arrives at 



a 2 

exp(W)W = -j exp(g 2 /4). (31) 

Since the Lambert function is defined as 

exp(W(s))W(s) = s, (32) 

one concludes that the location of the horizons as a function of q = \e\/M, is given by the 
real branches of the Lambert functions 

Aq 2 

X+ ~ ~4W(0,-^exp(g 2 /4))-g 2 ' (33) 

and 

4o 2 

x- = • (34) 

4W(-l,-^exp(g 2 /4)) - q 2 

The functions W(0,s) and W(—l,s) are the only real branches of the Lambert function 
with the branch point at s — — 1/e, where e is the base of natural logarithms. The horizons 
r + and r_ for 

\e\/M = 2W 1/2 (l/e), (35) 



merge at 



XM = 4Iy (f> , ( 36) 
extr l + W(l/e)' 1 ; 



where W(s) is a principal branch of the Lambert function W(0, s). Numerically one has 

x extr = 0.871, (37) 

and 



\e\ 
M 



= 1.056. (38) 



Inspection of (33) and (34) shows an interesting feature of the ABG geometry: the black hole 
solution exists for q greater than the analogous ratio of the parameters of the RN solution. 

The location of r + and r_ as a function of q for the charged black holes of both types are 
displayed in Fig. 1; its numerical values for some characteristic values of \e\/M are presented 
in Table II. Inspection of the figure shows that locations of the event horizons of the RN 
and ABG solutions are almost indistinguishable for, approximately, \e\/M <■ 0.7, whereas 
the differences between the inner horizons are more prominent. The latter differences are 
irrelevant here as in our analyses we shall confine ourselves to the static region exterior to 
the event horizon. Generally, for a given q, r + of the RN black hole is always greater than 
r + of the ABG geometry. 
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We remark here that the global structure of the ABG spacetime is similar to that of 
RN, with one notable distinction. Simple analysis shows that the curvature invariants of the 
curvature tensor, R, Ricci 2 , and Riem 2 are regular as r — > 0, and, moreover, other differ- 
ential invariants of the Riemann tensor and its covariant derivatives also exhibit regularity 
there. One concludes therefore that the ABG geometry for q < q extr represents the regular 
black hole solution. While this property of the ABG solution is not surprising it should be 
remembered that earlier efforts have been in unsuccessful in this regard. 

By means of the Wick rotation one obtains the Euclidean version of (25) with (28), which 
has no conical singularity provided the time coordinate is periodic with the period given by 



Ph 



Air lim 

r— »r+ 



U{r)V~ 




(39) 



Making use of elementary properties of the hyperbolic functions one has 



Ph 



An 



r±M 



AM J 



(40) 



We recall also that analogous period of the Euclideanized RN geometry is given by 

„2\-l 
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2e l 
T 



r 



(41) 



In the limit e —>■ both (40) and (41) tend to the Schwarzschild value 8nM whereas in 
the extremality limit Ph tends to infinity. As the Hawking temperature is proportional 
to the inverse of the period Ph one concludes that in the extremality limit the Hawking 
temperature of the ABG black hole vanishes. Moreover, closer analysis indicates that for 
a given e and M the ABG black hole is hotter than its RN counterpart characterized by 
the same values of the parameters. Of course, as expected, for small electric charges both 
temperatures are practically indistinguishable. 



IV. RENORMALIZED STRESS-ENERGY TENSOR IN THE SPACETIME OF 
ELECTRICALLY CHARGED BLACK HOLE 

In this chapter the method described in Sec. II is used to construct the renormalized 
stress-energy tensor of the quantized massive scalar fields in the ABG and extremal ABG 
spacetimes in the region exterior to the event horizon. As there are important simplifications 
for £ = 1/6 we shall consider only the conformal coupling. 

The analogous tensor in the RN geometry has been evaluated in Ref.[8] by means of 
the sixth-order WKB approximation of the solution to the scalar field equation and the 
summation thus obtained mode functions by means of the Abel-Plana formula. This result 
has been rederived and extended to the case of other spins, using simplified version of Eqs. 
(19-24) valid in the spacetimes with vanishing curvature scalar |fl6[ . 

Calculating the components of the Riemann tensor, its contractions and required covari- 
ant derivatives, inserting the results into (9-18), performing appropriate simplifications, and 
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finally constructing the renormalized stress-energy tensor, after rather lengthy calculations 
one has 



_L(T/) ren = - {576 7 M 7 r 6 (-626 7 M + 285 r) 



- 576 S e 2 M r (170 7 M — 669 7 Mr + 270 r z 



+ 485 e 4 M 4 r 4 [l86 (-47 + 181 (3) 7 2 M 2 

+ (3967 - 38132 (3 + 34165 /3 2 ) Mr + 8520 (3 r 2 

+ 125e 8 M 2 r 2 [2 7 2 M 2 (3431 (3 2 + 47307 /3 3 - 2976 - 32264 f3) 

+ Mr (2824 + 52712 (3 - 58927 (3 2 - 79068 /3 3 + 82459 (3 4 ) 
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965 e 6 M 3 r 3 
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= 4032 7 M 6 (22 7 M — 15 r) r 6 



+5 2 e 12 7 [(270 - 558/3 - 1500 1 + 2060 /3 3 + 1350 (3 4 - 1622 /3 5 ) M 
-15(9 -50 /3 2 + 45 /3 4 ) r 

+ 125 e 2 r { - 48 M 4 r 4 (478 7 2 M 2 - 391 7 M r + 54 r 2 ) 

+4e 2 M 3 r 3 [6 (101 + 873 (3) 7 2 M 2 

-(707 + 3148 (3 - 3855 /3 2 ) Mr + 472 (3r 2 

+e 6 Mr [2 (-544 - 1304 f3 + 957 (3 2 + 2041 /3 3 ) 7 2 M 2 

+r (392 + 408 (3 — 1717 /3 2 — 612 (3 3 + 1529 (3 4 ) M + 40 (3r (2 — 3 (3 2 ) 



-8e 4 M 2 r 2 



6(-73 + 35/3 + 222/3 2 ) 7 2 M 2 



-r (170 - 387/3 - 510/3 2 + 727 /3 3 ) M + 11 (3 (3 2 - l)r s 



+e 5 



2 (18 + 407 /3 + 441 /3 2 - 48 7 (3 3 - 551 /3 4 ) 7 2 M 2 



(42) 
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and 



7a 



+r (4 -287(3- 30 (3 2 + 676 (3 3 + 30 (3 4 - 393 /3 5 ) M 
+2 (2 - 15/5 2 + 15/3 4 ) r 2 ]}, 



(T e e ) ren = 576 7 M 7 (734 7 M- 315 r) r 6 



where 



+48e 2 M 4 r 4 
+ 1976 /3r 2 



+5e 2 [ - 576 M 5 r 5 (2202 7 2 M 2 - 1313 7 Mr + 162 r 2 ) 

6 (1399 + 3777/3) 7 2 M 2 + r (-3353 - 10914/5 + 14267 /3 2 ) M 
- 12e 6 M 2 r 2 [2 (-202 - 1198/3 + 749 /3 2 + 2161 /5 3 ) 7 2 M 2 
+r (154 + 1184 /3 - 2227 (3 2 - 1776 (3 s + 2665 /3 4 ) M + 92 (3 (-2 + 3 (3 2 ) r 5 
-48e 4 M 3 r 3 [6 (-91 + 649 (3 + 1042 /3 2 ) 7 2 M 2 



+ r (873 - 2308 (3 - 2619 /5 2 + 4054 (3 3 ) M + 206 



+12e 8 Mr 



3/3 2 )r s 



2 (80 - 519/3 - 1169 /3 2 + 1089 /3 3 + 1755 /3 4 ) 7 2 M 2 
+r(-176 + 639 (3 + 1320 /3 2 - 2716 (3 3 - 1320 /3 4 + 2253 (3 5 ) M 
+20 (2 - 15/3 2 + 15/3 4 ) r 2 " 

-e 10 [2 (129 + U1(3- 1590 /3 2 - 2750 (3 s + 2 1 05 (3 4 + 3001 /3 5 ) 7 2 M 2 
+3r(-43 - 272 (3 + 781 /3 2 + 960 (3 s - 2005 /3 4 - 720 (3 5 + 1299 /5 6 ) M 
+ 4 /3 (17 - 60 /3 2 + 45 /3 4 ) r 2 ]}, 

or 1 = 967r 2 m 2 x 60480M 5 r 15 , 



(43) 



(44) 



(45) 



/3 = tanh 



2Mr' 



(46) 



and 



7 = 1-/5, 



(47) 



(48) 



Obtained tensor is, as expected, covariantly conserved, and as could be easily verified in the 
limit e = 0, it reduces to the stress-energy tensor constructed in the Schwarzschild spacetime 
by Frolov and Zel'nikov. 

Equations (8-18) may be employed also in the RN geometry. Since the scalar curvature 
is zero there, both T^ ab and T^ ab vanishes, and the resulting tensor exhibits simple linear 
dependence on the coupling constant. Indeed, repeating the calculations for the line element 
(26), one has [8,16] 
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<«„ = (": - U l -)n, do) 



where 



Cf = - - ( 1248 e 6 - 810 r V + 855 M 2 r 4 + 202 r 2 e 4 

* 30240 7i 2 m 2 r 12 V 

- 1878 M 3 r 3 + 1152 MrV + 2307 MVe 2 - 3084 rMe 4 ) , (50) 

CI = \ - - (444 e 6 - 1488 Mr 3 e 2 + 162 r 4 e 2 + 842 r 2 e 4 - 1932 rMe 4 

30240 7T 2 m 2 r 12 V 

+ 315 M 2 r 4 + 2127 M 2 r 2 e 2 - 462 M 3 r 3 ) , (51) 



and 

= " ono , n I o 19 (3044 r 2 e 4 - 2202 M 3 r 3 - 10356 rMe 4 
30240 7?^ m z r LZ v 

+ 3066 e 6 - 4884 r 3 Me 2 + 9909 r 2 M 2 e 2 + 945 M 2 r 4 + 486 r 4 e 2 ) . (52) 

Since we are interested in the conformally coupled massive scalar fields, the exact form of 
the D b a tensor is irrelevant. 

Now, we shall address the question of how the differences between the geometry of the 
black hole spacetimes constructed within the framework of the Einstein-Maxwell theory on 
the one hand and the nonlinear electrodynamics coupled to the General Relativity on the 
other, are reflected in the overall behavior of the components of the stress-energy tensors. 
To answer this, let us analyze numerically (T^) ren in both cases. The results of our calcu- 
lations are presented graphically in Figs. 2-10. The plots of the time, radial, and angular 
components of the stress-energy tensor of the quantized massive scalar field as a functions 
of the (rescaled) radial coordinate in the spacetimes of ABG and RN black holes for three 
exemplar values of the ratio |e|/M = 0.1, 0.5and0.95, are supplemented by similar plots 
drawn for the extremal black holes. Inspection of the figures indicates that there are strik- 
ing qualitative similarities between the RN and ABG solutions for a given q. Moreover, for 
small values of the ratio the curves are practically undistinguishable from each other, and, as 
expected, noteworthy differences occur only for the black holes at and near the extremality 
limit. Since at large distances the line element (28) approaches that of the RN, the most 
interesting region is the neighborhood of the black hole event horizon. From Eq.(19) we 
know that the renormalized stress-energy tensor depends on the coupling constant £ in a 
complicated way, and, therefore, one should not expect that such similarities occur also in 
a general case. 

Specifically, the dependence of (T/) ren constructed in the spacetime of the nonlinear black 
hole and the RN geometry on r/r + for |e|/M = 0.1,0.5,0.95, are shown in Figs. 2 and 3, 
respectively. In Fig. 4 similar curves are drawn for the extremal black holes. In the most 
interesting region, i.e., in the vicinity of the event horizon, the energy density, p = — (T}) ren , 
is negative and decreases with increasing of |e|/M. 

As is seen in Figs. 5 -7, the radial component of the stress-energy tensor is everywhere 
positive, and the horizon values of the radial pressure, p r = (T^) ren , increases with increasing 
|e|/M. 
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Of all components of the renormalized stress-energy tensor, the most complicated be- 
havior exhibits the angular pressure pg = (Tjj) ren (Figs 8 - 10). Indeed, for the ABG black 
hole the angular pressure is positive on the event horizon for \e\/M < 0.937 and negative 
for larger values of the ratio. Moreover, for q < 0.903, (T{j) ren has a maximum at r = r + , 
whereas for larger values the angular pressure has its maximum away from the event horizon. 
Similarly, for the RN black hole pg is positive for q < 0.927 and it has its maximum away 
from the event horizon for q > 0.864. 

It could be checked by a direct calculation that 



lim 



<r;\ 



2M 



r / rein j 



tanh ■ 



2Mr 



(53) 



remains finite at the event horizon. We observe that since the DeWitt-Schwinger approxi- 
mation is local and the geometry at the event horizon is regular, one expects that that the 
renormalized stress-energy tensor is also regular there. 

It should be stressed once again that for arbitrary curvature coupling one has to incor- 
porate also the terms T^ ab , T^ ab , and T^ ab , that may considerably modify the results. 
Moreover, inspection of the Table I shows that for the neutral massive spinor and vector 
fields in the ABG spacetime one has to use the full system (9 - 18) while in the geometry of 
the RN black hole, the terms (9) and (11) do not contribute to the final result. 



V. CONCLUDING REMARKS 

In this paper we have constructed the renormalized stress-energy tensor of the massive 
conformally coupled scalar fields in the spacetime of the electrically charged black hole, 
being the solution of the coupled Einstein equation and the equation of nonlinear electrody- 
namics. A regular solution of this type has been recently given by Ayon-Beato and Garcia. 
The method employed here is based on the observation that the first order effective action 
could be expressed in terms of the traced coincidence limit of the coefficient a 3 . The gen- 
eral {T b ) ren , which has been obtained by functional differentiation of the effective action 
with respect to a metric tensor, has been applied in the spacetime of the nonlinear black 
hole. Since the Reissner-Nordstrom line element and ABG solution are practically indistin- 
guishable from each other for small values of \e\/M, one expects that this similarity should 
be reflected in the behavior of the renormalized stress-energy tensor. Explicit calculations 
confirm this hypothesis and show that important differences between appropriate tensors, 
(Ta)ren, evaluated in the spacetime of the RN black hole and that of ABG occur, as expected, 
near the extremality limit. For small q constructed tensors are practically indistinguishable. 
Moreover, analyses of the Hawking temperatures indicate that for a given mass and electric 
charge, the ABG black hole is hotter than its RN black hole counterpart. Since notable 
differences appear for temperatures close to zero one can ascribe this to the different ways 
of approaching the extremality limits. 

Apart from obvious extensions of our results to the massive scalar fields with arbitrary 
curvature coupling and to fields of higher spins, let us mention an interesting and important 
direction for future work. It is a problem of the back reaction of the quantized fields upon 
spacetime geometry of the ABG black hole, which may be studied perturbatively by means 
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of the semiclassical Einstein field equations with a source term given by the renormalized 
stress-energy tensor of the quantized massive field and the classical stress-energy tensor of 
the background nonlinear electromagnetic field. To guarantee the renormalizabilty at that 
level, the semiclassical equations should contain higher derivative geometric terms. It is 
especially important in view of the recent claim that the semiclassical zero temperature RN 
black holes do not exists [14]. 

It should be stressed that the DeWitt-Schwinger expansion is local, and, therefore, does 
not describe particle creation which is a nonperturbative and nonlocal phenomenon. The 
method also breaks down in strong or rapidly varying gravitational fields, and, moreover, 
the massless limit leads to the nonphysical divergences. However, it is expected that for 
sufficiently massive scalar field the DeWitt-Schwinger approximation provides a good ap- 
proximation of the exact renormalized stress-energy tensor. 
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r/M 




FIG. 1. The location of r + (upper branches) and r_ (lower branches) of the RN and ABG 
geometries as a function of e/M. The curve representing horizons of the ABG black hole is shifted 
to the right with respect to the one wich has been determined in the RN spacetime. 
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FIG. 2. This graph shows the radial dependence of the rescaled component {Tf) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field in the ABG geometry. From top to bottom the curves are for \e\/M = 0.95, 0.5, and 0.1. 
In each case (Tf} ren has its positive maximum at r = r + and attains negative minimum away from 
the event horizon. 




FIG. 3. This graph shows the radial dependence of the rescaled component (T/) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field in the RN spacetime. From top to bottom the curves are for \e\/M = 0.95, 0.5, and 0.1. 
In each case (T/) ren has its positive maximum at r = r + and attains negative minimun away from 
the event horizon. 
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FIG. 4. This graph shows the radial dependence of the rescaled component {Tf) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field. Top to bottom the curves are respectively for the extremal ABG geometry and the 
extremal RN black hole 



X<T T T > Ren 




FIG. 5. This graph shows the radial dependence of the rescaled component (T r r ) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally cou- 
pled scalar field in the Ayon-Beato and Gracia geometry. From top to bottom the curves are for 
\e\/M = 0.95, 0.5, and 0.1. In each case (T^) ren has its positive maximum at r = r+ and decreases 
monotonically with r. 
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FIG. 6. This graph shows the radial dependence of the rescaled component (T r r ) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field in the RN spacetime. From top to bottom the curves are for \e\/M = 0.95, 0.5, and 0.1. 
In each case {Tf) ren has its positive maximum at r = r + and decreases monotonically with r. 




FIG. 7. This graph shows the radial dependence of the rescaled component (T^) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field. Top to bottom the curves are respectively for the extremal ABG geometry and the 
extremal RN black hole 
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FIG. 8. This graph shows the radial dependence of the rescaled component (Tg) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field in the ABG geometry. From top to bottom the curves are for \e\/M = 0.1, 0.5, and 0.95. 
For |e|/M < 0.903, {T^) ren has its positive maximum at r = r+. For larger values of the ratio it 
approaches its maximum away from the event horizon. For \e\/M = 0.937 the angular pressure 
vanishes on the event horizon. 
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FIG. 9. This graph shows the radial dependence of the rescaled component (T[j) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field in the RN spacetime. From top to bottom the curves are for \e\/M = 0.1, 0.5, and 0.95. 
For \e\/M < 0.864, (T[j} ren has its positive maximum at r = r + . For larger values of the ratio it 
approaches its maximum away from the event horizon. For |e|/M = 0.927 the angular pressure 
vanishes at r + . 




FIG. 10. This graph shows the radial dependence of the rescaled component (Tg) ren 
[A = 90(8M) 4 m 2 7r 2 ] of the renormalized stress-energy tensor of the massive conformally coupled 
scalar field. Top to bottom the curves are respectively for the extremal ABG geometry and the 
extremal RN black hole. 



TABLE I. The coefficients a\ for the massive scalar, spinor, and vector field 



a{ 

0) 

0) 

0) 
a\ 

to 



2 S 5 S ^ 56 



5^ 
1 
140 



a 



a 



a 



a 



30 ( 6 £J 



Qn 



Q' 



10 



945 
2 
315 
1 

1260 
17 
7560 

1_ 

270 



1/2 



3 
140 
J_ 
14 
1 
432 
_ J_ 
90 
7 
720 
25 
376 
47 
630 
19 
630 
29 
3780 
_ J_ 
54 



27 
280 

_9_ 

28 
_5_ 
72 

31 

60 
_ J_ 
10 
_52 
63 
19 
105 

61 

140 
67 

2520 

J_ 

18 
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TABLE II. Location of r + and r_ of ABG and RN black holes for exemplar values of \e\/M. 



\e\/M r^/M (ABG) r + /M (ABG) r_/M (RN) r + /M (RN) 

0l6 0.001 1.995 0.005 1.995 

0.50 0.060 1.866 0.134 1.866 

0.95 0.422 1.356 0.688 1.312 

extremal 0.871 0.871 1 1 
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